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A  Propagator  Expansion  Method  For  Solving 
Linearized  Plasma  Kinetic  Equations 

With  Collisions 

1.  INTRODUCTION 

In  1946,  Landau1  gave  a  method  for  solving  linearized  plasma  kinetic 

problems  where  discrete  particle  interactions  were  neglected.  Subsequent 
2  -4 

studies  of  collisional  plasmas  have  employed  methods  of  solution  tailored 
to  a  specific  form  for  the  collision  operator.  In  this  paper,  we  give  a  general 
expansion  method  for  solving  linearized  plasma  kinetic  problems  when  colli¬ 
sions  are  included.  The  method  can  be  applied  to  a  wide  class  of  collision 

operators  and  has,  for  example,  produced  closed-form  results  for  the  colli- 

5  6 

sional  dielectric  function  for  the  Balescu  -Lenard  collision  operator. 

The  essence  of  the  idea  presented  in  this  paper  is  the  derivation  of  a 
collisional  propagator  in  terms  of  the  collisionless  propagator  for  the  Vlasov 
equation  and  the  linearized  collision  operator,  and  the  representation  and  use 

(Received  for  publication  21  June  1984) 
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of  the  collisional  propagator  either  as  a  series  solution  in  powers  of  the  collision 
frequency  v  or  as  an  asymptotic  expansion  as  the  collision  frequency  tends  to 


2.  KINETIC  EQUATIONS 

The  kinetic  equations  in  the  electrostatic  approximation  for  a  collisional 
plasma  are: 

la  _a  a  I  _  _  . 

I  at  +  v  ’^+L^re  E<r,  t>  +  a0(v)J  •  j  f(r,  v,  t)  =  £  (f)  ,  (1) 

~  •  E(r",  t)  =  4 jr  Jn^qj  +  q e  J d3v  ffr’.v’.tlj  ,  (2) 

where  E  is  the  self-consistent,  electrostatic  field;  ag  is  the  force -per-unit  mass 
due  to  spatially  uniform,  stationary  fields;  a, id  eC  (f)  is  the  nonlinear  collision 
operator  that  acts  on  the  velocity-space  coordinates  of  the  electron  distribution 
function  f.  The  ions  are  assumed  to  provide  a  fixed  neutralizing  background  with 
charge  density  n^qj  while  the  other  quantities  have  their  usual  meanings.  We 
seek  a  perturbation  solution  of  the  form 

f(r,v,t|  =  f0(v)  +  fjfr'.v',  t)  ,  (3) 

E(r,  t)  =  Ei<r,  t)  ,  (4) 


where  fo  is  the  zero-order  solution  and  fi  and  E^  are  perturbed  quantities.  The 
first-order  equations  using  abbreviated  notation  are: 

IS  .  9  .  -  <L_  .1/1  .1/  # . « 


(ft  +  v  •  ^  +  ^  <t)  =  S^t)  , 

9_  .  JL.  .  .  4  Jrqe  f  j3v  ' 

dr  dr  J 


dr  9r  J 

sv(t)  =  5_  ^(t)  •  L.  fo  .  (7) 

iti  dr"  dv 

e 

where  the  subscript  1  on  the  perturbed  quantities  has  been  deleted  and  the  super¬ 
script  v  is  used  to  denote  collisional  dependence,  (t)  is  the  electrostatic  poten¬ 
tial,  and  L1'  denotes  the  linearized  collision  operator.  The  operator  equation  for 
the  collisional  propagator  operator  Ul  (t-t'),  is 


where  U^(t-t')  is  given  by 


U^t-t')  =  fd3r'fd3v'Gv{r-r\v,v',t-t')  , 

and  where  Gv(?-r', v, v', the  Green's  function  for  the  problem,  is  transla- 
tionally  invariant  in  space  and  time.  Here,  Uv(t)  is  subject  to  the  causality  con 
dition  that  Uv(r)=0  forr  <0  and  the  initial  condition  that  Ul,(0)=l.  In  terms  of 
U^it-t')  the  solution  for  fv(t)  is 


fv(t)  =  Uv(t)  fv 


dt'  Uv(t-t')  SVW)  , 


(H 


where  fv(0)  denotes  the  initial  value  of  fv(t).  It  is  convenient  to  transform  Eqs 
(6)  through  (10)  to  k-v-w  space  using 


fV(t)  =  TiTjU‘-  ^.v.d]  =  f^u  )  ,  (1 

fWu  )  =  T£t  ^(r.v.t)  =  T^t  fv(t)  ,  (i: 


r  ~  iiot)  , 


(i: 


r-*'/A  l 


dt  exp(-ik  •  r  +  iwt) 


(I‘ 


where  k  is  real,  Im  ui  b>0,  and  the  contour  C  runs  from  -  oo+ib  to  +»+ib. 
Here  b  is  chosen  sufficiently  large  enough  to  ensure  that  the  t-integration  con- 
verges.  The  equations  in  fT-v-w  space  are: 


tv(u)  =  U^jV^fo)  + 


(i 


7.  Titchmarsh,  E.C.  (1948)  Introduction  to  the  Theory  of  Fourier  Integrals, 
Oxford  University  Press,  London,  Chapter  1. 


(16) 


■  4~i?  'V'01  • 


„  / 4)fqe  \  c  •?  -  a 

ev(u)  =  1  -  i  - -  / dJv  U^w)  k  •  —  fn  , 

\mek2/‘/  0 

where  fvj*(0)  is  the  spatial  transform  of  fl(0).  The  operator  Vv(u)  is 
Uv(u)  =  y*d3v'  Gv(k,7,v',u)  , 
and  the  Green's  function  in  k-v-u  space  is 


G^fk.v,  v'.io)  =  „  G^R.v,  v',t  )  . 

JK,  T 


(17) 


(18) 


(19) 


3.  EQUATION  FOR  THE  COLLISIONAL  PROPAGATOR 

In  order  to  derive  an  equation  for  U^fu),  we  consider  a  collisionless  propa¬ 
gator  operator,  U(t-t'),  which  obeys  the  operator  equation 


5=  ♦  r0'5y  )«<*-•'> -0 


(2  0) 


and  is  subject  to  the  same  causality  and  initial  conditions  as  those  imposed  on 
Ul  (r  ).  We  return  to  Eq.  (8)  and  write  an  operator  equation  for  Ul  (t-t1)  in  terms 
of  U(t-t')  and  Lv : 


u'(t-t')  =  U(t-t')  +  fX  dt"U(t-t")  Lr  U''(t"-t’> 

J  t' 


(21) 


where  we  have  used  causality  and  the  initial  conditions.  Transforming  Eq.  (21) 
to  k-v-w  space,  we  obtain 


Ul’(u)  =  U(w)  +  U(w)  Lr  Ur(u)  . 
Solving  Eq.  (22)  formally,  we  obtain 
~  ,,  -1  ~ 

U'  (w)  =  (  1  -  U(u)  L1  ]  U(u)  . 


(22) 


(23) 


Equations  (15)  through  (17)  and  Eq.  (23)  give  a  formal  solution  to  the  collisional, 
initial -value  problem. 


10 
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4.  ASYMPTOTIC  EXPANSIONS  AS  v  -  0 


In  order  to  pursue  a  solution  by  the  expansion  method,  we  introduce  a 
collision  frequency  v  and  write  Lv  =  v  L.  In  most  cases,  v  will  be  small  com¬ 
pared  to  other  frequencies  of  interest.  For  a  wide  class  of  collision  operators 
including  integral,  differential,  and  integral -differential,  and  for  some  con¬ 
tinuity  and  integrability  conditions  on  the  functions  that  appear  in  the  problem, 
it  can  be  shown  that  the  iteration  of  Eq.  (22)  yields  an  asymptotic  expansion  for 
Vl'(u)  as  t  -0; 

Ul'(u)  ~  U(w)  +  v  U(u>)  L  U((j)  +  . . . 

+  vn[U(u)  L]n  U(u)  +  .  .  .  .  (24) 

By  this  we  mean  that  if  g(u)  =  g(k,  v,  u)  and  fy(u)  =  Uv(u)  g(u>),  then  the  action  of 

the  iterated  U*'(u),  given  by  Eq.  (24),  on  g(u)  generates  an  asymptotic  expansion 

for  f^fu)  as  v  —  0.  Similar  results  can  be  obtained  for  the  iterated  ev  (w  ).  If  the 

iterated  fr(u)  is  a  divergent  asymptotic  expansion,  then  Eq.  (24)  may  not  be 

truncated  and  used  to  approximate  the  initial-value  problem  as  t-»  do  because 

time  secularities  will  appear  in  the  solution  when  transformed  to  r’-v’-t  space. 

However,  even  if  the  iterated  e^(u)  is  also  divergent  but  r/up«l,  where  Wp  is 

the  plasma  frequency,  then  e1'  (u)  may  be  truncated  at  first  order  in  v  and  set 

equal  to  zero  to  obtain  an  approximate  solution  for  the  roots  of  the  dispersion 

relation.  This  provides  a  solution  for  the  normal  mode  and  stability  problems 

accurate  to  order  v/u>  .  The  expression  for  the  collisional  dielectric  function 
P  -* 

truncated  to  the  first  order  in  v ,  denoted  by  t  (k,  w),  is 

2 

U) 

e(1)(k,w)  =  1  -  i  ^  {T0(k,u>)  +  v  T^k.u)}  .  (25) 

T  (k,  u)  =  f l3v  [  U(w)  L] "  U(u  )  k  •  —  fofv")  ,  (26) 

J  9v 


where  has  been  normalized  to  one. 


5.  UNIFORMLY  CONVERGENT  EXPANSIONS  IN  POWERS  OF  v 


For  the  special  case  when  U(u>)  is  the  free -particle,  collisionless  propagator 
and  L  is  an  integral  operator,  Eq.  (22)  is  an  integral  equation  for  Uv(u).  In  our 
study  of  problems  of  this  type,  we  have  found  that  if  the  functions  that  appear  in 
the  problem  are  continuous,  integrable,  and  obey  some  reasonable  uniform 
bounds,  then  the  iteration  procedure  yields  the  Neumann  series  for  the  associated 
integral  equation  and 

00 

Ur(u)  =  U(u)  +  £  rn[U(w)L]n  U(u)  ,  (27) 

n=l 

is  the  power-series  solution  for  U^u),  uniformly  convergent  in  the  neighborhood 
of  t/=0.  In  this  case  the  iterated  eI/(ui)  is: 


2  oo 
Up 


w)  =  1  -  i  —  'y  ^  i/nTn(kf  w)  , 
^  n=0 


X  k  ■  —  fo(v”)  , 
dv 


(28) 


(29) 


where  fQ  has  been  normalized  to  one  and  the  expression  for  the  free -par  tide, 
collisionless  propagator  has  been  used: 

U(u)  =  /dV  a  (v -v1)  ,  (30) 

*'  u-k-v 

For  problems  of  this  type,  it  also  follows  that  Eq.  (28)  is  the  power-series 
solution  for  e''(w),  uniformly  convergent  in  the  neighborhood  of  v  =  0.  If  Eq.  (27) 
can  be  summed  in  closed  form  and  analytically  continued  for  Im  w<b,  then  the 
initial-value  problem  is  solved. 


6.  SYMMETRY  PROPERTIES 


It  can  be  shown  that  G1  (k,v,v',u  )*  =  G!'(-k,v,v',  -w^Iwith  similar  formulae 
— ►  — ♦ 

for  the  other  k-v-u  quantities  discussed  in  Sections  2  and  3,  and  that  the  roots  of 
the  collisional  dispersion  relation,  e"(u)=0,  obey  the  relation  ui*’(k)  =  -u/'i-k)*. 
The  same  symmetry  properties  hold  for  each  of  the  iterated  k-v-ui  quantities 
discussed  in  Sections  4  and  5  when  truncated  to  the  Nth  order  in  v. 


7.  SOME  APPLICATIONS 

For  the  free -particle,  collisionless  propagator,  three  types  of  linearized 
collision  operators  have  been  studied:  integral,  differential,  and  integral-differ¬ 
ential. 

As  an  example  of  the  integral-type  collision  operator,  we  studied  the  simple 
BGK“  collision  operator,  given  by 

L1'  =  r  [-1  +  f0tv)Jd\'  ]  ,  (31) 

where  fo  has  been  normalized  to  one.  We  found  that  the  iterated  U'’(w)  and  the 
iterated  g'  (cj)  are  uniformly  convergent  power-series  solutions  in  r  that  agree 
with  expansions  of  the  exact  results. 

The  linear  simple  Fokker -Planck  collision  operator  is  a  differential  operator 
and  has  been  studied  by  Lenard  and  Bernstein.  It  is 

L1  =  v  ■  (  v  +  Vp  ■“)  ,  (32) 

9v  '  1  dv  ’ 

where  v-p  is  the  electron  thermal  speed.  We  found  that  the  iterated  U  (u)  is  a 
divergent,  asymptotic  expansion  as  i •  — 0  in  the  sense  described  above,  yet  the 
iterated  el (u)  for  fg  Maxwellian  is  a  uniformly  convergent  power-series  solution 
in  r  and  agrees  with  the  expansion  of  the  closed-form  result  for  e1  (u).  For 
example,  the  solution  to  the  dispersion  relation  in  the  long  wavelength  limit  to 

order  v  yields  the  total  damping  rate  which  is  the  sum  of  the  Landau 

L  K  ,, 

(collisionless)  part  y  jj*  and  the  collisional  part  y  where 

7r-  =  -  i  [  1/2  +  2(k/ko>2  +  ...  ]  •  (33) 

Here  kD  is  the  Debye  wavenumber.  Note  that  these  results  are  obtained  without 
solving  a  differential  equation  involving  Lr. 


The  linearized  Balescu'  -  Lenard  collision  operator  is  an  integral-differ¬ 
ential  operator.  We  have  found  that  the  iterated  Uv(o>)  is  a  divergent,  asymptotic 
expansion  as  v  —  0  and  we  have  obtained  closed-form  results  for  the  iterated 
tv(u)  to  the  first  order  in  v.  Using  these  results,  we  have  found  the  solution  for 
the  collisional  dispersion  relation  in  the  long  wavelength  limit  for  the  linearized 
Balescu  -  Lenard  collision  operator  and  for  fg  Maxwellian.  We  obtain  the  total 
damping  rate  r1  £■  which  is  the  sum  of  the  Landau  (collisionless)  part  and 
the  collisional  part  71  ,  where 

2 

\ 

D 


Ik 


4 


2  ^ei  f1  _2(  kD)  +  •"] 


(34) 


1/2 


no 


1/2  3/2 


in  A 


(35) 


I  /a  -* 

Here  rei  =  2  i  ee,  k  =  |k  I,  kD  is  the  Debye  wave  number,  n0  is  the  unperturbed 
density,  Te  is  the  unperturbed  temperature,  and  In  A  is  the  Coulomb  logarithm. 
Equation  (34)  is  new  and  was  obtained  in  collaboration  with  B.  Basu.  A  detailed 
discussion  of  its  derivation  is  quite  lengthy  and  will  be  published  elsewhere. 


a  DISCUSSION 

Equations  (21)  through  (23)  are  the  main  results  of  this  paper.  They  give 
equations  for  the  collisional  propagator  in  terms  of  the  collisionless  propagator 
and  the  linearized  collision  operator.  An  iterative  solution  of  Eq.  (22)  or  (23), 
given  by  Eq,  (24),  which  is  what  we  call  the  expansion  method  of  solution,  when 
substituted  into  Eqs.  (15)  through  (17)  provides  a  solution  of  the  linearized  plasma 
kinetic  equations  for  a  given  collision  operator.  This  method  of  solution  as 
applied  to  the  collisional  plasma  problem  discussed  here  is  new  and  is  particularly 
useful  when  a  direct  solution  of  the  problem  cannot  be  found.  For  example,  the 
problem  of  finding  the  direct  solution  of  the  linearized  plasma  kinetic  equations 
with  the  Balescu -Lenard  collision  term  and  then  using  it  to  find  the  collisional 
damping  of  the  electrostatic  plasma  waves  appears  to  be  intractable.  But  by 
applying  the  expansion  method,  we  have  been  able  to  obtain  the  collisional  damping  rate 
for  fg  Maxwellian  to  the  first  order  in  the  plasma  parameter,  which  is  given  by 
Eqs.  (34)  and  (35). 
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